Nonuniversality of the intrinsic inverse spin-Hall effect in diffusive systems 
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We studied the electric current induced in a two-dimensional electron gas by the spin current, in 
the presence of Rashba and cubic Dresselhaus spin-orbit interactions. We found out that the factor 
relating the electric and spin currents is not universal, but rather depends on the origin of the spin 
current. Drastic distinction has been found between two cases: the spin current created by diffusion 
of an inhomogeneous spin density, and the pure homogeneous spin current. We found out that in 
the former case the ISHE electric current is finite, while it turns to zero in the latter case, if the 
spin-orbit coupling is represented by Rashba interaction. 
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I. INTRODUCTION 

The spin-Hall effect (SHE) and the inverse spin-Hall ef- 
fect (ISHE) can be observed in two and three-dimensional 
electron systems with a strong enough spin-orbit inter- 
action (SOI)ii^. Via this interaction the electric current 
induces a flux of spin polarization flowing in the perpen- 
dicular direction and vice versa. These effects take place 
in metals and semiconductors, where the spin-orbit inter- 
action arises from impurity scattering, or band structure 
effects. Nowadays they are intensively studying theoret- 
ically (for a review see H) and experimentallyii^. These 
phenomena establish an important connection between 
spin and charge degrees of freedom that can be employed 
in spintronic applications. 

Here we will focus on ISHE. This effect is driven by the 
spin current which can be produced by different ways. It 
can be created by diffusion of an inhomogeneous spin po- 
larization, or can be induced directly by a motive force 
of various nature®'''. In experimental studies the former 
method was used in Refs. 5,8, while the latter was em- 
ployed in Refs. 1 91)10, From the theoretical point of view 
there are two quite distinct mechanisms of ISHE, de- 
pending on the extrinsic or intrinsic nature of SOI in an 
electron system. The extrinsic effect is promoted by the 
spin-orbit scattering of electrons from impurities^. The 
intrinsic effect is associated with the spin-orbit splitting 
of electron energy bands. This effect has been studied 
in Refllll together with the extrinsic mechanism. A sur- 
prising result of this study is that the finite inverse SHE 
takes place even in the case of a pure intrinsic Rashba^J. 
SOI, while the direct effect has been shown to vanish in 
the considered case of a diffusive systeni^'^. A reasonable 
explanation is that the Onsager relation between direct 
(SHE) and reciprocal (ISHE) effects should not be sat- 
isfied, because the spin-current is not conserving. This 
argument also means that for ISHE effect the coefficient 
in the local linear dependence Ic = Clg of the charge 
current density Ic from the spin current density Is can 
depend on the source that originally excites Is- In this 
sense ISHE is not universal. At the same time, SHE is a 



universal effect, because the coefficient relating Ig to Ic 
does not depend on how the electric current is produced. 
It can be created, for example, by electron diffusion, as 
well as by an external electric field. The result will be the 
same. It follows from the gage invariance of the electro- 
magnetic field. Formally, one obtains the same spin cur- 
rent, independent on whether it is induced by the scalar 
electric potential, or time dependent vector-potential. 

In order to demonstrate the non-universality of ISHE 
we will consider two kinds of spin-current sources. In the 
first case, an inhomogeneous spin polarization parallel to 
the z-axis creates the spin flux due to spin diffusion. In 
the second case, the spin current is driven by a spatially 
uniform "electric" fields, such, that the fields acting on up 
and down spins have opposite signs. The latter situation 
corresponds to spin current generation mechanisms sug- 
gested in Refs. Our goal is to show that the factors 
C are different in these two situations. Since our analy- 
sis has shown that in the case of the Rashba spin-orbit 
interaction C = for the source of the second kind, we 
will consider the cubic Dresselhaus interaction, as well, 
and demonstrate that the Onsager relation holds in this 
case. 

The outline of the paper is as follows. In Sec. II the 
linear response equations relating the spin and charge 
currents to the auxiliary fields will be written for a disor- 
dered two-dimensional degenerate (2DEG) electron gas. 
From this pair of equations the auxiliary fields can be 
excluded and linear relations between the electric and 
spin currents can be established. In Sec. HI this the- 
ory will be applied to the cases with Rashba (Subsection 
A) and Dresselhaus (Subsection B) spin-orbit couplings. 
The discussion of results will be presented in Sec. IV. 



II. LINEAR RESPONSE THEORY 

The Hamiltonian of the electron system has the form 

H = Ho + V, (1) 

where Hq is the unperturbed Hamiltonian of 2DEG, that 
includes the electron's spin-orbit coupling and their scat- 
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tering on randomly distributed spin-independent elastic 
scatterers. The spin-orbit coupling has the general form 

Hso = hk • (T, (2) 

where the effective magnetic field hk is a function of the 
electron momentum k and cr = (ax-, o'y^az) is the vector 
of Pauli matrices. In general, hk can be generated by 
the bulk-inversion asymmetry in the bulk and structure- 
inversion asymmetry in a quantum well (QW).^"* The 
perturbation term V — Vi + V2 represents interactions 
of electrons with the auxiliary fields. We will consider 
two types of fields. The first one is a slowly varying in 
time nonuniform Zeeman field B which is directed per- 
pendicular to 2DEG (z- direction). The corresponding 
interaction Hamiltonian is 

Vi^a.B. (3) 

Another interaction is 

V2 = cT,k ■ A, (4) 



This Hamiltonian contains the uniform spin-dependent 
field (TzA, where A slowly varies in time. Such a field 
induces the spin current by driving in opposite direc- 
tions electrons having opposite spins. It can be created, 
for example, by applying a time-dependent strain to a 
nonccntrosymmetric semiconductor. Indeed, as knowi*^^ 
the strain field u^z gives rise to the spin-orbit interac- 
tion aazUxzkx- Hence, in this case = au^z- Other 
mechanisms^ii of creating homogeneous spin currents can 
also be presented in a form of a spin dependent vector 
potential that is able to drive spins. 



Within the linear-response theory the spin P and 
charge I"^ currents of noninteracting electrons can be 
written in terms of retarded GJ^j^, (w) and advanced 
k' i^) single-particle Green's functions. Due to im- 
purity scattering these functions are nondiagonal with 
respect to the wavevectors k and k' . The linear-response 
expressions for the currents, as functions of the frequency 

and wavevector q, at 51 ^> are given by 
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where the spin-current and charge-current operators have 
the conventional fornJ^ j^/j, = (1/2)[ts/cj + , with v = 
k/m* -I- 9(hk • cr)/d'k and r, = az,Tc = e, np{uj) is the 
Fermi distribution. In the following the low-temperature 
case will be assumed, so that nF(j-o+^) — np{uj) — fl6(u)). 
The angular brackets denote averaging over disorder. 
This averaging will be performed within the semiclassic 
approximation, according to the standard procedure^'', 
where we will neglect the weak-localization corrections. 
We will assume that the spatial variations of the exter- 
nal field are slow within the electron mean free path Z, so 
that Iq ^ 1. This case corresponds to the diffusion ap- 
proximation, implying the expansion of Eq. ([5]) in powers 
of q. Also the SOI field will be assumed weak enough, so 
that hkp ^ 1/t, where r is the mean electron scattering 
time. 



III. INVERSE SPIN HALL EFFECT 

A. Rashba SOI 

Let us first consider ISHE in the case of Rashba spin- 
orbit interaction, where the spin-orbit field is linear in 
k and has the form hk = h^ = ak x z. If the auxil- 



iary field is Vi, given by Eq. ([3]), it creates a nonequi- 
librium and nonuniform in space spin polarization Sz- 
This distribution of electron spins relaxes to the uni- 
form state via diffusion, that is accompanied by a pure 
spin current. When V{fl,q) in Eq. ([5]) is represented by 
Vi(J7,q), the last term in this expression vanishes. Also, 
the terms containing the products CC and CG" can 
be shown to vanish, at least up to linear in q terms. Since 
in the following the higher-order terms starting from q^ 
will be ignored, only the products of the form G^G"' will 
be retained in Eq. ([5]). We assume that B in Eq. 
varies in x-direction, so that /"^ is expected to flow in 
y-direction. In Fig. 1 the Feynman diagrams contribut- 
ing to Eq. ([S]), where V = Vi, are shown. The upper 
(lower) arms in the diagrams denote the impurity aver- 
aged functions G^°'\uj) = (cj — iJk — hk • cr ± iT)~^ and 
the dashed lines depict the random scattering potential 
correlator (|C/kP). For simplicity this correlator will be 
assumed short-ranged, i.e. independent on fc, so that 
F = TriVp-dC/kp) = t:Nf\U\^ = l/2r is simply a con- 
stant. The multiple scattering blocks in diagrams shown 
at Fig. [IJb) and (d) represent processes where the initial 
electron spin density Sz evolves in the diffusion process to 
Si. Since this process is accompanied by the spin preces- 
sion due to Rashba SOI, i can be either z or x, as follows 
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from the spin diffusion equation'^^ for the spin polariza- 
tion varying in space along the x-coordinate. In general 
such a diffusion-precession dynamics is represented by 
the diffusion propagator Dij{c[). In the matrix form it 
can be represented as 



[(l-|C/P*/2)- 



(6) 



(7) 



takes a simple form in the case when q <C am*, that 
is when spatial variations of the Zeeman field are slower 
than spin-density variations caused by spin precession in 
the SOI field. In this case it follows from Eq. ^ that 
Dxz — \U\'^'^xzDxxDzz/'^- A straightforward calculation 
using Eqs. (jl]),© and ([T2|) gives for the first term in 
the brackets of Eq. ((TT|) the expression iTTNpqvpDzz/^'^ , 
that is twice larger and has opposite sign with respect to 
the second term. Finally, from Eqs. (ITO| and pT|) the 
charge current becomes 



Using the above definition, the contribution of all four 
types of diagrams in Fig. [T] can be written as 



ly i^B [Ky,D,, + KyxDxz + a—^Dxz , , 



where 



k 



TO* ^ 



(8) 



(9) 



It is easy to see that the diagonal components of D are 
finite at g — > 0, while the nondiagonal ones vanish as 
the first power of q. Therefore, in the leading approxi- 
mation the correlator K in the second term of Eq. 
must be calculated at q = 0. Up to the small semiclas- 
sic corrections of the order of {akp / Ep)^ this correla- 
tor is Kyx = —2TTaNp/r, and the last two terms can- 
cel each other. At the same time, it is easy to see that 
Kyy is at g = 0. Therefore, we did not include the 
corresponding term KyyDyz into Eq. ([S]). Further, as 
follows from Eq. ©, the correlator Ky^ is proportional 
to hk X hk+q. Therefore, it turns to at 5 = 0. In 
the leading approximation one finds from Eqs.(|H]El) that 
Ky:, = -i7rqa^NFkj,/2m*r^ and 
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qBD, 



(10) 



Our goal is to get an expression of the charge current 
through the spin current /|. Therefore, the next step is 
to calculate the spin current induced by the perturbation 
Baz- This current can be written in the form 

/| = ^^B{R-,,Dx.+ 

Dzzy^^Tr[azGl^jLo)azGU^)]] , (11) 



where 



K 



E 

k 



-Tr[a,Gl+^iL,)akGt{io)] . (12) 



The second term in the brackets of Eq. (|lip is equal 
to —iTTNpqvpDzz/2r'^ . This term represents the diffu- 
sion spin-current. In its turn the first term is associ- 
ated with spin precession caused by the Rashba field. It 
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(13) 



This result coincides with Ref. |ll|, taking into account 
that 2r = l/r and that the definition of /J in RefEll 
differs by the factor 1/2. 

The next example is the charge current induced in the 
y-direction by the external perturbation given by Eq. ^ , 
where A is parallel to the cc-axis. In this case the last 
term in Eq. ([S]) turns to zero, along with the terms con- 
taining the products CC and C^G"^. Further, a simple 
inspection of diagram (a) in Fig. 1 shows that it is zero 
at g = 0. The contribution of other diagrams to ly can 
be expressed as 



9 ^y^+ dk- 



(14) 



where the first term corresponds to Fig. 1 (b), while the 
second one is given by Fig. 1 (c) and (d). Since q — 0, 
only diagonal components of D enters in Eq. I^^. Also, 
at g = only i — x must be retained in the sum. As a 
result, after calculation of Kyx, one can see that that the 
sum in brackets turns into zero, up to the small semiclas- 
sic corrections of the order of {akp / Ep)^ ■ Therefore, 
within the semiclassic approximation the homogeneous 
pure spin current can not induce ISHE. At the same time 
the spin current created by this source is finite and is 
given by the Drude formula 



2r 



(15) 



This expression does not depend on the spin-orbit cou- 
pling. The latter enters as a small correction ^ ^feT^. 

Our calculations in this subsection show that ISHE is 
not universal. The induced by this effect electric current 
is finite, or zero, depending on whether the spin-current 
is produced by diffusion of an inhomogeneous spin polar- 
ization, or it is a pure uniform spin flux created by an 
external force of the form Eq. ^ . The driving force of 
this sort could be taken into account within the formal- 
ism employed in Ref. [l^ We, however, can not directly 
see if their expressions for spin and charge currents give, 
as we expect, vanishing ISHE, because these equations 
are presented in a rather general form. 
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the form 
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V 



I — A 
2tt 
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(17) 




(d) 




where 



and 



(18) 



FIG. 1: The Feynman diagrams for the charge current gen- 
erated by the intrinsic spin-Hall effect. The auxiliary field V 
can be either V\, or V2, where V\ and V2 are defined by Eq. Q 
and Eq. Q, respectively 



Dresselhaus SOI 



It is easy to see that the first term in Eq. p7|) turns to 
zero, similar to Eq. in the Rashba case. However, 
other two terms are finite, while they vanish for Rashba 
SOI, as well as for any other SOI which depends linearly 
on k. Taking SOI in the form of Eq. (|16p . from definitions 
(UHl), (Ull), (HH) and dS])-© one obtains at g = 



Although at y = V2 and for SOI given by the Rasha 
interaction the electric current is zero, we do not expect 
that the same takes place for a cubic in k Dresselhaus 
SOI. The reason is that the spin-Hall effect does not van- 
ish in the latter casC'^^. The Dresselhaus SOI field in a 
quantum well grown along the [001] direction is given 
by2i 



hl = pk,{kl 
— [iky{n1 



(16) 



where denotes the operator — (d/dz) averaged over 
the lowest subband wave function. Since h]^ is a nonlin- 
ear function of k, Vk/i^ entering into Eq. (fT4| is not a 
constant. Therefore Eq. (|14p has to be modified. Denot- 
ing by a bar the average Vk/i^ over the Fermi surface, 
the modified expression for the current can be written in 



Nf-, 







r2 \ dky ^ 



dky 



dky 



D:, 



vv 



T 



frf'l 



(20) 



Since only Ry^ is finite in Eq. PTI) . one has to calculate 
Qy. From Eq. (|18l) it can be expressed as 



Collecting all together one obtains from Eq. (|17p 



J 



Np 

ly = ^^^^^ 



2 "k"'^ 



dky 



■h^kx 



(22) 



This electric current can now be expressed through the 
spin current. The latter is induced by the time-dependent 
"vector potential" A in (jj]) and is given by Denoting 
by Q the expression in the square brackets of Eq. (|22]) . 



one obtains 



2eQ 
Tm*Vp 



(23) 
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Taking into account that Q oc /i^ it easy to see that the 
charge-to-spin current ratio is of the same order of magni- 
tude as in the case considered in Subsection A, Eq. ([T3l) . 
provided that the Rashba and Dresselhaus interactions 
are comparable in their strengths. One more useful re- 
lation can be obtained by using the expression for the 
spin-Hall conductivity derived in Ref. [20||2^ This con- 
ductivity can be written as ush = eNpQ/T^ . Expressing 
Q in Eq. p3)) through asH, and writing the electric con- 
ductivity in the form of the Einstein relation a = 2NpD 
we find 



On the other hand, the spin current induced by the spin- 
Hall effect is given by /| = asnE, where E is the electric 
field in the y-direction. Writing it as £' = /I* /cr we arrive 
to /J — (TsHly/cr. This equation, together with Eq. (p4)) 
establish Onsager relations between spin and charge cur- 
rents. 



IV. CONCLUSIONS 

Our analysis shows that the proportionality coefficient 
in the linear relation between the electric and spin cur- 
rent densities in the inverse spin-Hall effect depends on 
the origin of the spin current. Therefore, it is not pos- 
sible to introduce a universal parameter that determines 
a charge to spin current response. This non-universality 



is most clearly seen in the case of Rashba SOI, where a 
pure spin current produced by diffusion of an inhomo- 
geneous spin polarization gives rise to the finite electric 
current, while the latter is zero when the spin current is 
induced by a uniform in space force. In this situation, 
however, the ISHE produces a finite charge current, if 
SOI is represented by a nonlinear in k Dresselhaus SOI. 
It is important that in such a case the spin-Hall effect and 
ISHE obey the Onsager relation for coefficients relating 
the spin and charge currents. 

It should be noted that calculated above expressions 
for the spin and charge currents are related to local cur- 
rent densities, while experimentally measured are total 
electric currents, or electric potentials that are responses 
not to local spin currents, but rather to currents that 
are integrated over some distance (in 2D transport). For 
example, due to SOI the spin-current density created by 
spin diffusion oscillates and decay when the distance x 
from the spin-injection source is increasing. One has to 
integrate this current over x to obtain the total electric 
current induced by ISHE. Since the relation Eq. [T3] has 
the local form it will preserve after such an integration. 
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